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Unsteady aerodynamic loads were measured on an oscillating NACA 64A010 airfoil in the NASA Ames 11 by 
11 ft Transonic Wind Tunnel* Data are presented to show the effect of the unsteady shock-wave/boundary-layer 
interaction on the fundamental frequency lift, moment, and pressure distributions. The data show that weak 
shock waves induce an unsteady pressure distribution that can be predicted quite well, while stronger shock 
waves cause complex frequency-dependent distributions due to flow separation. An experimental test of the 
principles of linearity and superposition showed that they hold for weak shock waves while flows with stronger 
shock waves cannot be superimposed. 


Nomenclature 

C = chord of wing, 0.5 m 
C L = unsteady li ft coef ficient 
C Lm = static lift coefficient 

C L a = complex amplitude of unsteady lift coefficient per 

radian (nose up +) 

C Ma = complex amplitude of unsteady moment coefficient 
per radian (at leading edge, nose up + ) 

C p = static pressure coefficient, (P-PINF) /QINF 
C p a = complex amplitude of unsteady pressure coefficient 
per radian 

e =cos(a)f) +/sin(cof) 

/ - frequency, Hz 

k = reduced frequency, wC/2C/ 

M - freestream Mach number 

P - local static pressure, N/m 2 

PINF = freestream static pressure, N/m 2 

QINF = freestream dynamic pressure, N/m 2 

Re = chord Reynolds number 

U = freestream velocity, m/s 

a = complex amplitude of unsteady angle of attack 

a m - static angle of attack 

a> = radian frequency, a> = 2 rf 

Introduction 

U NSTEADY flows are one of the least understood 
problem areas in fluid mechanics. This is particularly 
true in the critical transonic flow regime where nonlinear and 
^dissipative effects tend to dominate the flow. Although recent 
advances in computational aerodynamics have made it 
possible to compute unsteady inviscid flows, unsteady flows 
with strong viscous interactions are not adequately modeled. 
Because of this lack of a computational model, the present 
understanding of transonic unsteady flows has been 
developed mostly through experimental studies. 

Transonic unsteady airfoil testing was initially motivated 
by the destructive effects of aileron buzz encountered on high- 
speed aircraft during World War II. Early work at the (then) 
NACA Ames Aeronautical Laboratory by Erickson and his 
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coworkers 1 * 2 on control surface oscillations culminated in the 
development and operation of a large test rig utilizing splitter 
plates in the (then) Ames 16-ft High-Speed Wind Tunnel. 
Unsteady loads on oscillating airfoils and control surfaces 
were measured in a series of tests over a five-year period. 3 7 
This early work at Ames was instrumental in developing 
miniature instrumentation and refining the newly acquired 
unsteady testing capability. Parallel work at the (then) NACA 
Langley Aeronautical Laboratory 8 * 10 clarified the effect of 
unsteady wall interference and wind-tunnel resonance. A 
summary of British and European work on oscillating airfoils 
in this period may be found in Refs. 11 and 12. All this early 
unsteady research was hampered at supercritical speeds by 
choking of the older unventilated wind tunnels. 

After almost 20 years of inactivity in two-dimensional 
transonic unsteady aerodynamics, there is currently an urgent 
technological need for a more complete understanding of 
these flows. The work of Tijdeman and his co-workers at the 
National Aerospace Laboratory (NLR), Amsterdam, The 
Netherlands, has motivated much of the current research. 
Using a carefully calibrated pneumatic-tubing system, they 
were able to extract first-harmonic unsteady aerodynamic 
pressures on models ranging from oscillating flaps to 
supercritical airfoils. 1315 Many of the new physical 
phenomena observed and reported in Refs. 13-15 are 
described in detail in Ref. 15. Recently, many other research 
groups have initiated oscillating airfoil experiments. Results 
of these investigations, mostly on smaller models, may be 
found in the proceedings of recent national and international 
conferences. 16 * 17 The current experimental investigations at 
Ames were designed to expand the existing unsteady 
aerodynamic test envelope to higher Reynolds numbers, to 
different modes of motion, and to more diverse mean flow 
conditions. 

The Ames oscillating airfoil experiment was run in the 1 1 by 
1 1 ft transonic wind tunnel. A 0.5 m chord x 1 .35 m span wing 
with an NACA 64A010 airfoil section was installed in the test 
section between two splitter plates as shown in Fig. 1. The 
model was driven into a pitching or plunging motion by a 
servocontrolled hydraulic actuation system. Motion was 
imparted to the wing via four push-pull rods imbedded in the 
splitter plates. The unsteady data were measured with an 
array of miniature in-situ piezoresistive pressure transducers. 
Steady data were measured with conventional pneumatic 
instrumentation. All the data were input directly into a 
minicomputer -based data system for subsequent online data 
reduction and analysis. More detailed descriptions of the test 
apparatus, wing model, instrumentation, and computer 
system may be found in Refs. 18-20. 
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HEIGHT/CHORD RATIO = 6.00 

Fig. 1 Oscillating airfoil installed in Ames 11 by 11 ft Transonic 
Wind Tunnel. 


NACA 64A010 



MACH NUMBER, M 

Fig. 2 Measured static-response curves for NACA 64A010 airfoil. 


The purpose of this paper is to discuss selected portions of 
the data that illustrate the complicated aerodynamic effects 
on a? oscillating airfoil at transonic speeds. An especially 
important component— one that is stressed in this paper— is 
the effect of the shock-wave/boundary-layer interaction on 
the unsteady loads. This effect is illustrated by examining in 
detail the first-harmonic unsteady pressure distributions at the 
two mean flow conditions shown in Fig. 2. Condition 1 
represents a transonic flow with a weak shock wave, while 
condition 2 is an example of a transonic flow with shock- 
induced separation (shock stall). For each of these conditions 
first-harmonic unsteady pressure distributions, generalized 
forces (lift and moment), and selected comparisons with 
theories are presented. These comparisons will show that 
weak shock waves induce an unsteady aerodynamic response 
that is amenable to present-day computational techniques. It 
is also shown that under shock stall conditions the unsteady 
pressures are not well predicted by an inviscid analysis and 
strong frequency-dependent viscous effects appear. The paper 
concludes with an experimental verification of the principles 
of linearity and superposition for these two conditions. These 
principles are routinely assumed in present-day practice and 



Fig. 3 Static pressure distribution in transonic flow (condition 1). 



Fig. 4 Effect of mean angle on static pressure distribution in 
transonic flow (condition 1). 


their limits of validity must be fully explored before 
meaningful transonic analyses can proceed. 

The data presented in this paper form a small part of the 
extensive data bank that has been assembled from this ex- 
periment. Some other data may be found in Ref. 19 and a 
First-harmonic analysis of the complete data set may be found 
in Ref. 21. Detailed information on the geometry of the 
NACA 64A010 test airfoil may be found in Refs. 22 and 23. 

Transonic Unseparated Flow 

This case is represented by condition 1 in Fig. 2 and lies just 
to the right of the boundary separating subcritical and 
supercritical flows. The supercritical flow region is not ex- 
tensive, and it is terminated by a weak normal shock. The 
flow is expected to be attached over the entire airfoil and, as 
will be shown, inviscid-flow modeling adequately predicts the 
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Fig. 5 Effect of frequency on unsteady pressure 
distributions. 


PITCHING 1 @ 0.25C 



subsonic theory. 


major features of the unsteady pressure distribution. The 
static pressure distribution is shown in Fig. 3, where a weak 
shock wave is shown at approximately midchord. Down- 
stream of the shock wave the mild adverse pressure gradient 
insures good pressure recovery to the trailing edge. 

The quasisteady response for the pitching mode is shown in 
Fig. 4. Since the motion is symmetric, only upper-surface 
pressures are shown. The magnitude of the unsteady pressure 
is approximately proportional to the distance between the 
a m = ±2 deg curves at each X/C station. Starting from the 
trailing edge, the response is minimal over the latter 25 
chord and increases rapidly in the vicinity of the shock wave. 
Upstream of the shock, the response is uniformly large. This 
behavior differs from that found in both subcritical ex- 
periments and linear theory where the unsteady amplitude 



x/c 


Fig. 7 Static pressure distribution compared with numerical 
solutions. 


increases smoothly from the trailing edge. 19 As will be shown, 
this flow condition is representative of Tijdeman’s Type II-l 
motion (Ref. 15, p. 72) since the weak shock wave disap- 
peared during a portion of the cycle. The shock-wave 
movement (where it exists) is consistent with inviscid-flow 
theory since it moves aft with increasing angle of attack. 

The measured unsteady pressure at the fundamental 
frequency is shown in Fig. 5 for a pitching mode. The real 
part is consistent with the quasisteady analysis given 
previously. The unsteady pressure is very small in the trailing- 
edge region, increases rapidly to a peak at the shock wave, 
and remains fairly level over the forward half of the airfoil. 
Its shape remains constant over the entire frequency range, 
but the amplitude decreases slightly with frequency. The 
imaginary part is out of phase with the motion upstream of 
the shock, but it quickly undergoes a 180 deg phase change at 
the shock. Unlike the real component, the imaginary com- 
ponent is locally antisymmetric with respect to the shock 
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position. Starting from a low level, the imaginary component 
increases with frequency to a level approaching that of the 
real component. The imaginary components have both 
positive and negative contributions to the loads. Also, unlike 
the real component, the imaginary component contributes 
significant postshock pressure disturbances. The unsteady 
pressures on the lower surface of this symmetric airfoil are 
always 1 80 deg out of phase with those on the upper surface. 

A comparison of the supercritical unsteady pressure 
distribution with linear subsonic theory at fc = 0.20 (f- 34 Hz) 
is shown in Fig. 6. The effect of the shock wave clearly in- 
validates the linear theory. Both the real and imaginary parts 
predicted by the theory are wrong. However, some interesting 
integrated results were obtained by using linear subsonic 
theory. When the areas under the real components were 
computed they were unexpectedly close. The real component 
of the lift and the center-of-pressure position are actually 
predicted quite well by using linear theory. (Notice that the 
“center-of-gravity” of the experimental curve is at 25 chord 
since the loading drops to zero downstream of the shock.) 
This result must be considered fortuitous and is probably 
accounted for by the favorable shock-wave position. The 
imaginary part is not well predicted, and even the real part 
must be quite sensitive to shock-wave position. 

A calculation was made by Magnus 24 for this flow con- 
dition by using the Euler equations of motion. Magnus 1 mean 
flow computation is shown in Fig. 7 along with the prediction 
based upon a steady transonic small-disturbance code 25 
(program TSFOIL). The predictions based upon Euler’s 
equation show good agreement with the data except for a 
slightly downstream shock-wave position. The small- 
disturbance theory predicts a slight overexpansion over the 
forward part of the airfoil relative to both the experiment and 
Euler’s equation. The adverse pressure gradients approaching 
the trailing edge are predicted quite well, an indication of fully 
attached flow over the entire airfoil. 

An unsteady pressure distribution from a harmonic analysis 
of Magnus’ data is shown in Fig. 8. The experimental data are 
the same as those presented in Fig. 6. Upstream of the shock 
wave, both the real and imaginary components are predicted 
quite well. The theoretical unsteady pressures caused by the 


PITCHING V @ 0 25C 
M = 0.8 
a m = 0° 

Re = 12 x 10 6 
k = 0.2 
f * 34 Hz 



Fig. 8 Upper-surface unsteady pressure distribution compared with 
numerical solution of Euler’s equation. 


shock motion are overpredicted, and the postshock response 
is also over predicted. This behavior is to be expected because 
the forward region of the airfoil is influenced mainly by the 
inviscid flow, while the shock dynamics and, to some extent, 
the postshock behavior must be influenced by the unsteady 
shock-wave/boundary-layer interaction. More extensive 
comparisons of these data with another solution of Euler’s 
equation that includes some viscous effects may be found in 
Ref. 26. Overall this result is very encouraging, and the 
outlook is that more efficient codes that include viscous ef- 
fects will be able to predict these flows quite well. 

The aerodynamic transfer function for this flow condition 
is shown in Fig. 9. The unsteady lift and moment is similar in 
many ways to that predicted by linear theory even though the 
flow is supercritical. The magnitude (fortuitously) agrees 
quite well with linear subsonic theory, but the phase lag of the 
supercritical data with respect to subsonic theory is quite 
large. This is caused by a combination of wake disturbances 
and some contributions from the shock wave. Scale effects do 
not seem to be important for this flow condition. For the 
reasons discussed previously, the predictions based upon the 
Euler equation overestimate the loads. 

Transonic Flows with Shock-Induced Separation 

The flow condition indicated by point 2 in Fig. 2 is 
representative of a transonic flow with a strong shock wave 
and subsequent flow separation. As will be shown, the 
oscillatory pressures exhibit many new characteristics that are 
not modeled by inviscid flow theory. The static pressure 
distribution is presented in Fig. 10. The shock wave is stronger 
and farther upstream than that shown in Fig. 3. The poor 
trailing-edge pressure recovery indicates probable flow 
separation. Interferograms obtained by Johnson and 
Bachalo 27 on a 15 cm chord model of the same airfoil clearly 


PITCHING ±l c AT 0.25 C 

M • 0 8 

~ Sr n° 


o R«- 12 X 10 6 


• EULER’S EQUATION 
LINEAR AERODYNAMICS 



coefficient, b) leading-edge moment coefficient (both per radian). 
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Fig. 12 Effect of mean angle on static pressure distribution for shock 
stall (condition 2). 


THEORY 



Fig. 11 Static pressure distribution compared with numerical 
solution of transonic small-disturbance equation. 


show an extensive separation region at these conditions. The 
importance of viscous interactions, even in the steady flow, is 
shown in Fig. 11 where steady transonic small-disturbance 
theory (code TSFOIL) is compared with the experimental 
data. An attempt was made to force better agreement by 
matching the theoretical lift to the experimental lift, but this 
attempt was unsuccessful. 

The quasisteady response for ±2 deg of pitch at a mean 
angle of attack of 4 deg is shown in Fig. 12. As mentioned 
previously, the magnitude of the unsteady response is 
proportional to the distance between the limiting angles. 
Starting with the lower surface, the response decreases with 
distance from the trailing edge, vanishes at about midchord, 


and is out of phase with the motion over the forward half of 
the airfoil, (the phrase “out of phase” means decreasing 
suction with increasing angle of attack.) Turning to the upper 
surface and starting from the trailing edge, the magnitude 
decreases to zero at about 63<Fo chord, is out of phase for 
about 20 °7o of the chord, vanishes again at 43 chord, and 
increases toward the leading edge. On the upper surface there 
are two points of vanishing response (A7C = 0.43 and 0.63). 
This curious behavior is caused by the unnatural shock-wave 
movement. Usually a shock wave moves downstream with 
increasing angle of attack. In this case, shock stall causes the 
shock wave to move upstream with increasing angle of attack. 
This effect has been observed for a long time in steady-state 
testing (Ref. 27, for example), but the shock stall has an 
extremely important effect on the oscillating airfoil. 

The unsteady pressures are presented in Fig. 13. The 
condition shown in Fig. 13, pitching ± 1 deg at 25 chord, is 
a somewhat less exaggerated version of the ± 2 deg excursions 
shown in the quasisteady response. At low frequencies the 
real part of the upper-surface response has both positive and 
negative lobes. The upper surface response is out of phase 
with the motion over the central portion of the airfoil. Both 
the theoretical and the experimental results shown in Figs. 6 
and 8 indicate no cases in which the in-phase response has this 
particular shape. The authors are not aware of any unsteady 
aerodynamic theories that can predict this behavior — what is 
needed is a mathematical model that includes strong unsteady 
viscous interactions. With increasing frequency, the real part 
of the response curve changes shape and, at high frequencies, 
again resembles inviscid flow. The lower-surface pressure also 
changes somewhat with increasing frequency, but not as 
dramatically as that of the upper surface. The imaginary 
component is also quite different from that of the previous 
case. A large bump first appears on the upper surface at 
moderate frequencies and persists to high frequencies. The 
lower-surface response is much less affected by increasing 
frequency. 

The predicted unsteady pressure distribution from linear 
subsonic theory is independent of mean angle of attack and is 
identical to that shown in Fig. 6; the linear theory is obviously 
inapplicable. 

The aerodynamic transfer function has been computed and 
is shown in Fig. 14. The experimental data are compared with 
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Fig. 13 Effect of frequency on unsteady pressure 
distributions. 
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Fig. 14 Aerodynamic transfer function for shock-stall condition: a) 
lift coefficient, b) leading-edge moment coefficient (both per radian). 


linear subsonic theory. Unlike Fig. 9, even the overall trends 
of the response are not well predicted by the linear theory. The 
magnitude of the loads peaks at about A: = 0.2 and is ac- 
companied by a rapid phase shift that is somewhat 
reminiscent of a one degree-of-freedom resonance. 

Linearity and Modal Superposition 

The principle of superposition is the foundation on which 
the theory of aeroelasticity is built. The underlying assump- 
tion that the generalized forces are linear functions of the 
generalized coordinates leads to the concept of superimposing 
two modes of motion to predict a third unknown mode. 



Fig. 15 Variation of unsteady lift coefficient magnitude with 
oscillatory angle of attack. 


The principle of superposition was tested by first measuring 
the loads for three modes, assuming linearity, and finally 
combining two of the measured modes to predict the third. 
The two superimposed modes were pitching at 50% chord and 
plunging. The predicted mode (also measured) was pitching at 
25% chord. The results are shown in Table 1. The un- 
separated flow tabulations show an error of less than 2% in 
magnitude and 2 deg in phase over the entire frequency range. 
This error is about the level of experimental accuracy, and 
may be taken as a good indication of the validity of modal 
superposition. The tabulation for shock-induced separation at 
two Reynolds numbers shows errors in predicted and 
measured loads that are greater than can be accounted for by 
experimental techniques alone. It is known that flutter 
boundaries, for example, can be very sensitive to small 
changes in the phase of the loading with respect to the motion. 
The data show that the predicted lift coefficient slightly leads 
the motion while the measured lift coefficient slightly lags the 
motion. It can be concluded that the superposition principle is 
violated in this case, but the question of whether or not 
superposition can be assumed for practical calculations must 
await definitive aeroelastic computations. 

The linearity of the dynamic loads was tested by varying the 
amplitude of the motion while keeping all the other 
parameters fixed. A representative example is presented in 
Fig. 15 for the magnitude of the lift coefficient. The attached 
flow, represented by a m =0 deg, follows linear theory closely, 
deviating only slightly at the higher oscillatory angles of 
attack. The separated flow, a m =4 deg, neither follows linear 
theory nor remains linear over an appreciable range of 
oscillatory angles of attack. This figure shows in another way 
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Table 1 Measured and predicted unsteady lift coefficient 
(NACA 64A010, M — 0.8, pitching at 25% chord) 


Measured Predicted* Error 

Rex 10 Mag Ph, deg Mag Ph, deg Mag, Vo APh, deg 


0.025 

12.6 

Unseparated flow, a m =0 deg 
0.171 - 8.8 0.169 - 9.0 

1.1 

0.2 

0.05 

12.6 

0.160 -16.1 0.159 -16.1 

0.6 

0.0 

0.10 

12.6 

0.134 -26.7 0.133 -28.3 

1.0 

1.6 

0.20 

12.6 

0.097 -27.2 0.095 -28.9 

1.9 

1.7 

0.20 

6.7 

0.091 -28.2 0.092 -29.5 

Shock induced separation, a m = 4 deg 

-0.7 

1.3 

0.20 

11.9 

0.166 - 0.9 0.153 4.6 

7.8 

-5.5 

0.20 

6.2 

0.126 - 4.1 0.132 3.3 

-5.4 

-7.4 


Linear superposition of measured pitching at 50^o chord + measured plunging modes. 


that separated flows are not necessarily superimposable. The 
deviations from linearity that are acceptable must again be 
found from practical computations and sensitivity analyses. 

Conclusions 

The data presented in this paper showed that transonic 
shock-wave/boundary-layer interactions have an important 
effect on unsteady pressure distributions and loads. For weak 
interactions, numerical computations are feasible and some 
progress has already been made in this direction. When the 
shock wave is strong enough to induce separation in the 
steady flow, the unsteady pressure distributions are strongly 
affected by the separation. Numerical modeling of steady- 
state separated flows is still not resolved, and the prediction of 
viscous, unsteady, transonic flows must await future 
developments. 

The question of modal superposition and linearity has also 
been investigated experimentally. The superposition principle 
has been verified for unseparated flows, but some deviations 
were measured for separated flows. Setting the limits to which 
these principles can be relaxed for practical transonic ap- 
plications must await detailed sensitivity calculations with 
realistic geometrical, elastic, and inertial parameters. 
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